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1 Time-Dependent Gross-Pitaevskii Equation 

The remarkable achievement of Bose-Einstein condensation in dilute trapped 
alkali-metal atomic gases flplpl has stimulated the (now successful) search for 
quantized vortices that are usually associated with external rotation H,H,pP,B| ■ 
An essential feature of these condensates is the order parameter, characterized 
by a complex macroscopic wave function ^(r, t). Theoretical descriptions of vor- 
tices in trapped low-temperature condensates have relied on the time-dependent 
Gross-Pitaevskii (GP) equation p jl0| , |ll| , which omits dissipation. For a trap 
rotating at angular velocity J? about z, it is a nonlinear Schrodinger equation 



ih^- = (T+V tr -f2L z +g\V\ 2 )V, (1) 



where T = — fi, 2 V 2 /2M is the kinetic-energy operator, V tI — \M^2,-tjj 2 x 2 is 
the harmonic trap potential, L z — xp y — yp x is the z component of angular 
momentum, and g = 4irh a/M characterizes the strength of the short-range 
interparticle potential (here, a is the positive s-wave scattering length for repul- 
sive two-body interactions; typically a is a few nm). Current experiments involve 
"dilute" systems, so that the zero-temperature condensate contains nearly all N 
particles, with f dV \&\ 2 w N. For a steady solution, &(r,t) = &(r)e~ illt/r \ 
where fi is the chemical potential. 

1.1 Equivalent Hydrodynamics of Compressible Isentropic Fluid 

If the condensate wave function is written as <?" = e |<^|, the real and imaginary 
parts of Eq. (fil) precisely reproduce the time-dependent irrotational hydrody- 
namics of a compressible isentropic fluid, written in terms of the number density 
n = |<f| 2 and the velocity potential $ = TiS/M, where v = W §Q|||. In 
principle, the dynamics of a curved vortex line in a rotating trap follows directly 
from general hydrodynamics, but the trap potential and resulting nonuniform 
density complicate the problem considerably. Thus it is preferable to start with 
simple situations. 
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1.2 Thomas- Fermi Limit for Large Condensates 

For a noninteracting condensate, the macroscopic wave function is the ground 
state of the harmonic oscillator, with spatial extent dj — (h/Mujj) 1 ^ 2 and 
j = x,y,z; typically dj is a few //m. The repulsive interactions act to expand 
the condensate, and the relevant dimensionless interaction parameter is Na/do, 
where do — (d x d y d z ) is a suitable geometric mean. Recent experiments fo- 
cus on the regime Na/do 3> lj when the quantum-mechanical energy associated 
with the density gradients is negligible. In this Thomas- Fermi (TF) limit fl3|| , 
the nonrotating GP equation has the solution g | &tf 1 2 + Kr ~ /^> with a parabolic 
density profile 

„„M- W <»)(l-|-^-£). (2) 

Here Rj — {2^,/MuJ 2 ) 1 / 2 ^$> dj fixes the condensate's dimensions and titf(0) = 
n/g is the central density. The (large) dimensionless ratios -Rq/^o = ^5Na/do 
and n/hujo = ^R^/d 2 , characterize the effect of the repulsive interactions, where 
ujq and R are appropriate geometric means. 

2 Energy of a Vortex in a Large Rotating Trap 

The time-dependent GP equation (fy) follows from a variational energy functional 

E{(2) = JdV<P*(T + V tI + y\>P\ 2 -(2L Z )V. (3) 

If W represents a straight singly quantized vortex in a rotating disk-shape trap 
(R± 3> Rz), the TF density yields the increase in energy AE(ro, Q) associated 
with the vortex at the transverse position r$ JL4|. For a nonrotating trap with 
J? = 0, this energy AE(tq, 0) decreases monotonically with lateral displacement. 
Energy conservation requires that the allowed motion is a precession at fixed trap 
potential (the trajectory is elliptical for an anisotropic trap). In the presence of 
weak dissipation, the vortex slowly spirals outward, lowering its energy. 

For arbitrary Q and small lateral displacements, AE(ro, SI) has the form 
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where £ = h/(2Mfi) 1 / 2 is the vortex-core radius and 2/R\ = 1/R 2 , + 1/R 2 y 
defines the mean transverse radius of the condensate. In the TF limit, £ is small, 
with £i? = dg, ensuring a clear separation of length scales £ <C d n «C Rq. 

• For small positive J?, AE{tq,Q) decreases with increasing lateral displace- 
ments so that a vortex will spiral outward in the presence of weak dissipation. 



Stability of a Vortex 



The curvature at ro — vanishes at fi m = ^(h/MR 2 l )ln.(R±/£), signaling 
the onset of metastability. Since AE(0, Q m ) > 0, this state is not truly 
stable. 

When fi > fi m , the trap center becomes a local minimum. For weak dissi- 
pation, a vortex that is slightly displaced from the center will spiral inward. 
When fi reaches fi c = |(^/Mi2j_)ln(i2x/0 = f^m, a central vortex be- 
comes stable because AE(0, fi c ) = 0. Hence fi c is the thermodynamic critical 
angular velocity for vortex creation. 

For fi > fi c , the form of AE(ro, fi) shows that a vortex at the outer edge 
remains metastable, but the barrier for entry into the condensate decreases. 
Spontaneous nucleation of a vortex may eventually occur through a surface 
instability |T^,|l6|,[l7l . The observed critical angular velocity for vortex cre- 
ation JJ3J7J is sa 70% higher than the predicted TF value fi c , in qualitative 
agreement with such a surface mechanism (but see Ref. |L8| and Sec. 4.2 for 
an alternative explanation). 



3 Small- Amplitude Excitation of a Vortex in a Rotating 
Trap 

A macroscopic Bose condensate acts like an external particle source, in the sense 
that the same physical excited state (here labeled by j) can be achieved ei- 
ther by adding or by subtracting one particle. The actual eigenstates involve 
"quasiparticle" operators a ■ and a,- that are linear combinations of the two 
quantum-mechanical states. For a given normal mode, the resulting pair of cou- 
pled complex amplitudes Uj(r) and Vj(r) obey the Bogoliubov equations |I(],|l9| 
that determine the corresponding eigenfrequency u)j. Imposing Bose-Einstein 
commutation relations [aj,aj,] = 5jk requires that these amplitudes obey the 
particular normalization J dV {\uj\ 2 — \vj\ 2 ) = 1, and the resulting quasiparticle 
Hamiltonian reduces to a set of uncoupled harmonic oscillators ^ ?iu;.,<xaj, 
summed over all modes with positive normalization. In the simplest case of a 
uniform condensate, the eigenfrequencies are all positive, which ensures that 
the system is stable because the energy then has a lower bound. For a uniform 
condensate moving with velocity i>, however, some of the eigenfrequencies can 
become negative as soon as \v\ exceeds the speed of sound, which corresponds 
to the instability associated with the Landau critical velocity [g(|. Physically, 
the system can spontaneously generate quasiparticles because the Hamiltonian 
is no longer bounded from below. 

These general ideas have direct relevance to the stability of a singly quan- 
tized vortex in a rotating trapped condensate. For simplicity, it is convenient to 
consider an axisymmetric condensate in equilibrium with a trap rotating at an 
angular velocity fi. In this case, states can be characterized by their azimuthal 
angular quantum numbers rrij, and the transformation to rotating coordinates 
H — > H — fiL z ensures that the eigenfrequencies ujj(fi) in the rotating frame 
are simply ujj(fi) = luj — mjfi, where u>j is the corresponding eigenfrequency in 
the nonrotating frame. 
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3.1 Stability of a Vortex 

The first numerical study of the Bogoliubov equations for a singly quantized vor- 
tex in a nonrotating trap found a single negative-frequency mode (the "anoma- 
lous mode") [ plp^ |, implying that the vortex in the condensate is unstable. 
Physically this instability arises because the condensate with a vortex has a 
higher energy than the vortex-free condensate. 

This situation is especially clear for a noninteracting Bose gas in an axisym- 
metric harmonic trap, when the first excited (vortex) state has an excitation en- 
ergy fku± and an angular momentum H. The transition back to the true ground 
state involves the (negative) frequency — u>± and the (negative) change in angu- 
lar quantum number — 1. More generally, the numerical analysis for small and 
medium interaction strength Na/d± |2l| found that the anomalous frequency u a 
remained negative (with m a = — 1 unchanged). In a frame rotating with angular 
velocity fi, the anomalous frequency becomes uj a (fi) = w a + fi. Since Lu a < 0, 
the eigenfrequency u} a (fi) in the rotating frame rises toward with increasing fi 
and vanishes at a critical rotation speed fi* = — lu u — \u a \- The Bogoliubov de- 
scription of small oscillations implies that a condensate with a singly quantized 
vortex is unstable for fi < fi* but becomes stable for fi > fi* . 

In contrast to the numerical study |gl) for small and medium values of the di- 
mensionless coupling parameter Na/do, a direct perturbation analysis is feasible 
for a large condensate (Na/do S> 1) containing an axisymmetric singly quan- 
tized vortex (the TF limit). Detailed study of the Bogoliubov equations Q 
yields the explicit expression for the anomalous frequency in the rotating frame 
uj a (fi) = -§ (h/MR 2 ± ) In (R±/0 + fi, which has the expected form. The result- 
ing critical angular velocity fi* = | (U/MRj_\ In (R±/£,) for the onset of local 
stability agrees with fi m inferred from Eq. (Q) for the onset of metastability. 

3.2 Splitting of Normal-Mode Frequencies Caused by a Vortex 

The ground-state condensate can sustain dynamical oscillations driven by the 
mean-field repulsive interaction (analogous to plasma oscillations in a charged 
medium) . These normal modes become particularly simple in the TF limit of a 



large condensate [23 , and experiments have confirmed the predictions for the 
lowest few modes in considerable detail [ |ll| . For an axisymmetric condensate, 
the normal modes can be classified by their azimuthal quantum number m, and 
modes with ±m are degenerate. 

When the condensate contains a vortex, the asymmetric circulating flow af- 
fects the preceding normal modes. In particular, the originally degenerate modes 
are split by the Doppler shift of the local frequency (analogous to the splitting 
of magnetic sublevels in the Zeeman effect). In the TF limit, this small frac- 
tional splitting of the degenerate modes is proportional to m|do/i?§ Q,g5|; it 
has served to detect the presence of a vortex and to infer the circulation and 



angular momentum [|7|,[26| . 
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4 Vortex Dynamics 

At zero temperature, the time-dependent GP equation (uh determines the dy- 
namics of the condensate in a rotating trap. A vortex line will move in response 
to the nonuniform trap potential, the external rotation, and its own local cur- 
vature. This problem is especially tractable in the TF limit, because the small 
vortex core radius £ permits a clear separation of length scales; the method of 
matched asymptotic expansions yields an explicit expression for the local velocity 



of a vortex line [|27 28 



4.1 Dynamics of Straight Vortex 

It is simplest to consider a straight vortex H], which applies to a disk-shape 
condensate with R± ;§> R z . Assume that the vortex is located near the trap 
center at a transverse position r (t). In this region, the trap potential does not 
change significantly on a scale of order £. The solution proceeds in two steps. 

(a) First, consider the region near the vortex core that is assumed to move 
with a transverse velocity Viz. Transform to a co- moving frame centered on 
the vortex core, where the trap potential exerts a force proportional to Vj_Vtr 
evaluated at r (t). The resulting solution includes both the detailed core struc- 
ture and the "asymptotic" region \r± — r± \ 3> £ 

(b) Second, consider the region far from the vortex, where the core can be 
treated as a distant singularity. The short-distance behavior of this solution 
includes the region £ <C \r± — r±o\. The two solutions must agree in the common 
region, which determines the translational velocity V of the vortex line. 

The details become intricate, but the final answer is elegant and physical: 



v= m 



AMix 



. R ± \ 2Mf2R 2 , 



W 3n 



Z X Vj.Vtr , (5) 



where R± for an asymmetric trap is defined below Eq. (|4|). This expression has 
several important aspects 

• The motion follows an equipotential line along the direction z x Vj_Vt r , 
conserving energy, as appropriate for the GP equation at zero temperature. 
For an asymmetric harmonic trap, the trajectory is elliptical. 

• For a nonrotating trap, the motion is counterclockwise, in the positive sense. 

• With increasing external rotation, the translational velocity V decreases and 
vanishes at the special value Q m = ^(h/MR 2 L ) ln(i?j_/£) discussed below 
Eq. (@). 

• For fl > firm the motion as seen in the rotating frame is clockwise. 

• A detailed analysis of the normalization of the Bogoliubov amplitudes shows 
that the positive-norm state has the frequency 

^i + ^y 
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This normal-mode frequency is negative (and hence locally unstable) for 
f2 < firm but it becomes locally stable for fl > £2 m , in agreement with the 
discussion below Eq. (Eh. 

4.2 Inclusion of Curvature 

Equation (|5|) for the local velocity of a straight vortex oriented along z can be 
generalized to include the possibility of a different orientation of the vector i 
locally tangent to the vortex core. In addition, local curvature k of the vortex 
line defines a plane that includes both t and the local normal n, inducing an 
additional translational velocity along the binormal vector b = t x n. A detailed 
analysis |29| yields the translational velocity of the element located at To 



T w s ft ftxVV tr (r ) , lt \, [ t \ 1 , k 2 \ , 2VKr(r )xfl 

V ( r o) = -7777 



2M \ 9 \^tf\ 2 "7 \^ R 2 ± ' 8 J ' VlV tT (r ) ' 

(7) 
where V^_ is the Laplacian in the plane perpendicular to fi. In the first term, 
|^tf| 2 vanishes near the condensate boundary; hence t x VVtr(^o) must also 
vanish there, implying that the vortex is locally perpendicular to the surface. 

Equation (|7|) allows a study of the dynamics of small-amplitude displacements 
of the vortex from the z axis, when x{z,t) and y(z,t) obey coupled equations. 
In the limit u z = 0, there is no confinement in the z direction, and the density 
is independent of z. The resulting two-dimensional dynamics exhibits helical 
solutions that are linear combinations of two plane standing waves. 

More generally, for oj z / 0, the density near the z axis has the TF parabolic 
form, and the solutions become more complicated. It is convenient to define the 
asymmetry parameters a = R^/R 2 and f3 — R 2 /R 2 Z1 where a > 1, (3 > 1 indicate 
a disk shape and a < 1, < 1 indicate a cigar shape. 

• For a nonrotating trap with the special asymmetry values a = 2/[n(n +1)] 
(here, n a positive integer), the effects of the nonuniform trap potential and 
the curvature just balance, and the condensate has stationary solutions with 
the vortex at rest in the xz plane. A disk-shape trap has no such states, and 
the first one occurs for the spherical trap with a = 1. The next such state 
occurs for a = i, when the condensate is significantly elongated. Similar 
considerations for f3 apply to stationary states in the yz plane. 

• For other values of a and (3, solutions necessarily involve motion of the vortex 
line relative to the stationary condensate. 

• Analytical solutions can be found for an extremely flat disk with a ^> 1 and 
/3 > 1, reproducing the frequency f2 m found in Eq. (^). 

• For small deformations of a vortex line in an axisymmetric trap with a = fi, 
a disk-shape or spherical condensate (a > 1) has only a single (unstable) 
precessing normal mode with a negative frequency uj a < 0. In this case, an 
external rotation J? > Q m = \u) a \ stabilizes the vortex. For these geome- 
tries, fl m is less than the thermodynamic critical value fl c . In a spherical 
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condensate, the one anomalous mode \uj a \ agrees with the observed vortex 
precession frequency seen in the JILA experiments [p||lq] . 
In contrast, an axisymmetric cigar-shape condensate has additional negative- 
frequency precessing modes, and fi m can exceed J? c for sufficiently elongated 
condensates; such behavior seems relevant for the ENS experiments ppjlq] , 
where fl m ~ 1.7 fl c is close to the observed rotation speed for creating the 
first vortex. 

For small axisymmetric deviations from a spherical trap, a straight vortex 
line can execute large-amplitude periodic trajectories. In this case, the vortex 
line becomes invisible when it tips away from the line of sight, and it then 
periodically returns to full visibility. Such revivals agree with preliminary 
observations at JILA |2ql. 
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